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ABSTRACT 

Studies of brown dwarf (BD) outflows provide information pertinent to questions 
on BD formation, as well as allowing outflow mechanisms to be investigated at the 
lowest masses. Here new observations of the bipolar outflow from the 24 Mjup 
BD, 2MASSJ12073347-3932540 are presented. The outflow was originally identified 
through the spectro-astrometric analysis of the [OI]/16300 emission line. Follow-up 
observations consisting of spectra and [SIX], R-band and I-band images were obtained. 
The new spectra confirm the original results and are used to constrain the outflow PA 
at ~ 65°. The [OI]/16300 emission line region is spatially resolved and the outflow 
is detected in the [SII] images. The detection is firstly in the form of an elongation 
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of the point spread function along the direction of the outflow PA. Four faint knot- 
like features (labelled A-D) are also observed to the south-west of 2MASSJ 12073347- 
3932540 along the same PA suggested by the spectra and the elongation in the PSF. 
Interestingly, D, the feature furthest from the source is bow-shaped with the apex 
pointing away from 2MASSJ12073347-3932540. A color-color analysis allows us 
to conclude that at least feature D is part of the outflow under investigation while A 
is likely a star or galaxy. Follow-up observations are needed to confirm the origin of 
B and C. This is a first for a BD, as BD optical outflows have to date only been de- 
tected using spectro-astrometry. This result also demonstrates for the first time that 
BD outflows can be coUimated and episodic. 

Subject headings: interstellar medium: jets and outflows - stars: pre-main- sequence 



Introduction 



The brown dwarf (BD) 2MASSJ1 2073347-3932540 (hereafter 2M1207A), with an estimated 



mass of 24 Mjup ( jMohanty et al.||2007| ), has garnered much attention in recent years. Not only 
does it exhibit young star-like properties, i.e. it is a strong accretor with a detected accretion disk 
( |Riaz et al.|2012b|at |Riaz & Gizis|2008[ [Mohanty et aL][2003| [2005] ) and bipolar outflow ( jWhelan 



et al. 2007), but it also has a planetary mass companion 2M1207B. 2M1207B was first detected 



by Chauvin et al. (2004) using direct imaging and it is one of only a small number of exoplanets 
to have been detected in this way. The 2M1207 system is a member of the TW Hya association 
(TWA). TWA is one of the closest young stellar associations to the Sun with an estimated distance 
of 47-67 pc and an age of ~ 10 Myr ( [Barrado Y Navascues||2006| ). Thus TWA is a key reference 
for the study of star, BD and planet formation and associated phenomena. It is the interesting 
properties of 2M1207 combined with its mass at the lower end of the BD mass spectrum and 
relatively close proximity, that have made this object the target of studies aimed at understanding 
BD and planet formation. 

In this paper new observations of the bipolar outflow driven by 2M1207A are presented. As 
studies of outflow properties can provide information on the central engine of the driving source, 
observations of BD outflows can contribute much to the debate on the origin and evolution of BDs. 
For example, the kinematics and position angle (PA) of the outflow can be used to infer the disk 
inclination and PA ( Whelan et al.||2009b ). Episodic or knotty outflows point to variable accretion 
activity and precession or wiggling in an outflow can be interpreted as evidence of an unresolved 
companion ( Whelan et al.||2010| ). In addition, measurements of the mass outflow to accretion rate 
{MoutlMacc) are relevant to jet launching models which in turn afi'ect models of disk evolution and 
thus planet formation ( Ferreira et al.|2006| Hartigan et al.|1995 1. One possible diff'erence between 
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BDs and low mass stars is that first results indicate that MoutlMacc may be higher for sub-stellar 
objects than low mass stars. For low mass stars it is currently estimated at ~ 10%, however, in BDs 
this ratio has been found to be significantly higher and, in the most extreme cases, the two rates are 
found to be comparable ^Bacciotti et al.|[20TT , |Whelan et al.||2009b ). This is significant as it may 
help to explain why BDs never accrete enough matter to become stars ( Machida et al.|20()9 ). Better 
observations of a larger number of BDs are needed to confirm these results, to properly understand 
BD accretion/outflow properties and thus to compare BD properties with low mass stars. 

BD outflows have mainly been studied through spectroscopy and especially the spectro- 
astrometric analysis of their forbidden emission line (FEL) regions. The first indications that 



BDs could drive outflows came from the detection of FELs in their optical spectra (Fernandez 
& Comeron 2001[ ). Spectro-astrometry (SA), which involves using Gaussian fitting to measure 
the centroid of an emission region as a function of wavelength/velocity, was used to confirm that 
these FEL regions were extended and originated in an outflow ( Whelan et al.|2005 ). Since the first 
detection, four further young BD candidates have been found to drive outflows using this method. 
Imaging studies of BD outflows have rarely been conducted and the mapping of the molecular out- 
flow from the BD ISO-Oph 102 with the sub-millimeter array (SMA), is the highest quality image 
in existence at present. The very low mass star (VLMS) Par-Lup3-4 with a mass just above the 
substellar limit ( Comeron et al.|[2003 ) has been studied with a combination of imaging and spec- 



troscopy, and a bipolar outflow was detected using both methods (Comeron & Fernandez 2011 



Fernandez & Comer6n||2005] ). Previous to this work it was the lowest mass protostar where an 
optical outflow had even been marginally detected using imaging. For Par-Lup3-4 the [SII] image 
is well extended along the outflow PA and a Herbig Haro (HH) object (HH 600) is also detected. 



Comeron & Fernandez (20 11 1 estimate a proper motion of 168 + 30 kms for HH 600. Wang & 



Henning (2006) carried out a survey for HH objects in the Chamaeleon I molecular cloud with the 
aim of searching for optical outflows from BDs. While 18 HH objects and ~ 40 BD candidates 
were identified, they found no persuasive evidence that any of the BDs were related to the HH 
objects. 

Optical spectra of 2M1207A obtained in May 2006 contained double peaked [OI]/16300 emis- 
sion which was tracing a bipolar outflow ( Whelan et al.||2007( ). This was an exciting discovery as 
it confirmed 2M1207A as the lowest mass galactic object known to drive an outflow. The observa- 
tions presented here can be divided into three separate results. Firstly, new optical spectra confirm 
the original detection and a comparison with the 2006 spectra allows an estimate of the outflow 
PA to be made. Secondly, we present a [SII] image in which the outflow is detected close to the 
source in the form of an elongation in the point spread function (PSF). Thirdly, a series of knot-like 
features are seen in the [SII] image along the PA suggested by the spectro-astrometric study. We 
explore the possibility that these features are shocks in the 2M1207A outflow. Below these new 
observations are described, their significance discussed and the best direction for future studies of 
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BD outflow activity considered. 



2. Observations and Data Analysis 



Both spectroscopic and imaging studies of the outflow driven by 2M1207A have been con- 
ducted with the European Southern Observatory's (ESO) Very Large Telescope (VLT), using the 
UV- Visual Echelle Spectrometer (UVES) and the FOcal Reducer/low-dispersion Spectrograph 1 
(FORS-1) dPekker et aL||2000t [Appenzeller et aL][2000| [Szeifert et al.|[Tj98| ). Follow-up obser- 
vations to the original UVES data were made with FORS 1 in January/March 2008 and UVES in 
August 2008. Starting first with the UVES 2008 spectroscopy, observations were carried out with 
the CD3 grating and a central wavelength of 6000 A, providing spectra with a wavelength range of 
4900-7100 A and spectral resolution R=40,000. The slit width and pixel scale were 1" and a.' 182, 
respectively, and the seeing was ~ 0'.'8. Two observations with exposure time totaling 2900 s were 
taken. The UVES 2006 data had the same instrumental set-up except that the central wavelength 
was set at 5800 A meaning that the wavelength range included the HyS line at 4861 A. Indeed, 
2M1207A has been observed numerous times with UVES and the spectra are available in the ESO 
data archive. However, observations were generally aimed at studying accretion and were not op- 
timised for outflow studies. This means that the exposure time was too low to detect the forbidden 
emission. 



The UVES 2006 observations, published in [Whelan et al.| ( |2007[ ), were made with the slit 
ahgned with a FA of 0°, i.e N-S. For the 2008 spectra the FA was set at the parallactic angle which 
varied between 75° and 85° over the course of the observation. By comparing the analysis of the 
2008 spectra with the 2006 spectra the FA of the outflow can be probed and, for the purpose of 
this we take the slit FA to be 80°. Figure [T] presents the difl"erent slit R^s used and the direction 
on the sky to which the offsets correspond to. The spectra were reduced using standard routines 
for bias subtraction, flat fielding and wavelength calibration, provided by the image reduction and 
analysis facility (IRAF). In order to analyse the spectra using SA, the 2D spectra corresponding to 
each individual order had to be extracted and wavelength calibrated using the Identify, Reidentify, 
Fitcoords and Transform routines. Continuum subtraction was carried out using the Continuum 



routine. Overall the process followed for the SA was the same as outlined in Whelan et al.|(|2007|) 



including the smoothing of the data using Gauss to increase the S/N ratio. Also see Whelan & 



Garcia (2008 1 for further details on the spectro-astrometric technique. 



Imaging with FORSl was carried out using the [SII]-l-62 narrowband filter and the R-BESS 
and I-BESS broadband filters in March 2008. All observations were made with the high resolution 
(HR) collimator and a binning of 2x2, giving a pixel scale of 0'.' 1 25 . Again standard IRAF routines 
were used for bias subtraction and flat-fielding. For the [SII] images a total of 2 x 30 mins expo- 
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sures were made. The exposures were taken consecutively, however the seeing deteriorated from 
~ 0'.'6 to (y.'9 between the two exposures. The total exposure time for the broadband images was 
15 mins. Although the [OI]/16300 line had previously been detected in the 2M1207A spectrum 
and shown to trace an outflow, for the purposes of imaging, [SIX] was preferred over [01] due to 
the relatively bright background emission in [01] from the atmosphere. To help with the analysis 
of the [SII] FORSl images, SUSI2 I-band images were obtained from the ESO archive and are 
further discussed in the Section 3. 

As dithering was not used to obtain the I-band images, fringing effects which could not be 
removed using the flat fields presented a problem to the analysis of faint emission in the vicinity of 
2M1207A. To remove the fringing, separate I-band observations were identified and downloaded 
from the ESO archive. The closest ones in time were taken on 1 1 February 2008, which is about one 
month before our observations of 2M1207, and had been originally obtained to image SN 2007 Y in 
the galaxy NGC 1 187. The galaxy has a relatively large angular size that fills a significant fraction 
of the field of view of FORSl. However, the area of the detector where 2M1207 was later imaged 
corresponded to a region in the distant outskirts of the galaxy with no obvious structure. Thanks 
to this, it was possible to construct a fringe pattern frame suitable for the region immediately 
surrounding 2M1207 by combining and median-filtering the eight available exposures of NGC 
1187, uncorrected for telescope offsets between the exposures. The fringe pattern obtained was 
visually confirmed to be very similar in appearance to that in the 2M1207 exposures. Next, an 
average sky background value near the position of 2M1207 in the detector was subtracted from 
both the 2M1207 and fringing pattern images. The background-subtracted fringing pattern image 
was then scaled by a numerical factor and subtracted from the background-subtracted image of 
2M1207. Based on visual inspection of the resulting image, this last part of the process was iterated 
with different scaling factors applied to the background-subtracted fringing image, minimizing the 
residual fringing pattern in the image of 2M1207. As dithering was used to obtain the SUSI2 
images of 2M1207 identified in the ESO archive, construction of a fringing pattern from separate 
images was not needed in this case. Instead, the frames containing the target were combined with 
median filtering to construct a sky frame, including fringing, to be subtracted from the individual 
images. The sky background remained very stable during all the SUSI2 observations, thus leaving 
no noticeable fringing residuals in the sky-subtracted images. 
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3. Results 



3.1. Spectroscopy and Spectro-astrometry 



In Figure [2] the ID extracted spectrum from the UVES 2006 data is shown. This spectrum 
was extracted using the ESO UVES pipeline and the purpose here is to show the full extent of the 
UVES spectra and to mark the strongest lines. The 2006 UVES spectrum is shown as the 2008 
spectrum does not include the region of the Hfi line. In this study SA is applied to the [OI]/16300, 
Ha and Hfi lines to search for evidence of the outflow and a comparison is made between the 2006 
and 2008 spectra. Beginning with the [OI]/16300 line, Figure [3] (a) compares the line profiles of the 
[OI]/16300 emission in the 2006 data and in the 2008 data. In both cases the line is double peaked 
with the blue and red-shifted peaks at ~ -8 kms~' and +4 kms"\ respectively. The small radial 
velocities of the peaks agree with the hypothesis that 2M1207A has a near-edge-on accretion disk 
and, therefore, an outflow close to the plane of the sky. The inclination of the 2M1207A disk has 
recently been estimated by Riaz et al. ( |2012b|a ) at 57° to 69°. The kinematics of the outflow and 
spectro-astrometric results would place the inclination at the upper end of this scale. In the 2006 
spectra, the emission peaks were found to be off'set in opposite directions by ~ 80 mas confirming 
the origin of the emission in an outflow. In Figure |3] (b) the position velocity diagram of the 
2008 [OI]/16300 emission region is shown. The two peaks are clearly seen and, interestingly, they 
are spatially ofl^set in opposite directions, verifying the conclusion of ,Whelan et al. ( 2007[ ) and 
revealing that the PA of the 2M1207A outflow lies closer to a PA of 90° than 0°. This is the first 
observation where a BD EEL region is spatially resolved and therefore where SA is not necessary 
to demonstrate origin in an outflow. In Figure |4] the PA of the outflow is recovered by plotting 
the offsets at 0° against the offsets at 80° as a function of systemic velocity. Taking into account 
the fact that the two slit PAs were not perpendicular the outflow PA is constrained at ~ 65° + 10°. 
Assuming that the outflow is perpendicular to the disk this gives a likely disk PA of ~ 155°. The 
PA of the companion is 125.6° + 0.7° ( [Mohanty et al.|2007p . The difference of ~ 30° between the 
disk PA (inferred from the outflow PA) and the PA of the companion could be accounted for if the 
companion and disk are not coplanar ( Riaz et aL]|2012b|a[ ). Alternatively, projection effects could 
explain the difference. The difference between the outflow axis projected on the sky and the PA of 
the planet will be 90° only at the time of the greatest elongations. Most of the rest of the time it 
will be near but not exactly 90° for a close to, but not exactly, edge-on orbital plane, depending on 
the position on the orbit. 

In addition to FELs, permitted emission lines are also strong tracers of outflow activity in BDs 
and low mass stars (Whelan et al. |2009a 2004). Indeed Ha and HjS emission is known to trace 
outflows and coUimated jets (Bacciotti et al. 2011[ ). Typically it is found that both outflow and 
accretion process contribute to the permitted line regions with the wing-emission often tracing the 
outflow component ( |Podio et al.]|2008| [Whelan et a'L]|2004| [Takami et aL]|2001D . SA is also used 
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here to search for any contribution to the Ha and Hfi lines from the 2M1207A outflow. [Whelan 



et al. ( 2007| ) analyzed the 2006 Ha emission region and no evidence of an outflow component 



was found. The Ha and Hfi lines of 2M1207A are, however, known to be highly variable (Stelzer 



et al.||2007 , Scholz et al.||2005p , and therefore one possibility that is important to consider, is that 



the 2006 spectra were taken during a period of high accretion activity, making it very difficult to 
detect any outflow component. If spectra are taken during a period of lower accretion activity the 
ratio between the emission from the outflow and the emission from the accretion could be more 
conducive to the detection of the outflow component. With this in mind, all the Ha emission from 
all the spectra taken in 2006 and 2008 were analyzed using SA. For the analysis of the [OI]/16300 
line the spectra were median combined to increase the S/N in the [OI]/16300 line region. However, 
the S/N in the region of the Ha line is high enough to analyse the spectra separately. Five individual 
spectra were obtained in 2006 with a further two in 2008. Figure [5] compares the Ha line profile 
shape and spectro-astrometric signature in all seven spectra. While the shape changes a lot (note 
the appearance of a small "bump" between the two main peaks), no spectro-astrometric signature 
is detected. Hence it can be argued that the majority of the emission is coming from close to the 
source, presumably from accretion, and that the Ha should provide a good estimate of the mass 
accretion rate of 2M1207A. The small "bump" is interesting as it was not previously detected in 
the variability studies of jStelzer et al. (2007) or Scholz et al. (2005|). Similarly we compare the 
line profile shape of the H/3 and spectra- astrometric signature of the HyS line. Again no spectro- 
astrometric signal is detected and only the line profiles are presented in Figure [6} 



3.2. [SII] Imaging: Elongated PSF 

The FORS-1 [SII] image with an inset of a zoom on 2M1207A is shown in Figure |7j Firstly 
we chose to look for evidence of emission from the outflow close to the BD by analyzing the image 
with the best seeing (~ 0'.'6). Interestingly, the PSF appears to be extended along the outflow PA 
suggested by the spectro-astrometric results. In Figure|9]the PSF of 2M1207A in the [SII] image is 
compared with the PSF of two nearby stars in [SII] and with 2M1207A in R. Each PSF was fitted 
with an ellipse and the results compared. The 2M1207A PSF is more elliptical than the others and 
is elongated along a PA given by the elliptical fit at 62° + 6°. This is consistent with the estimate of 
65° + 10° for the outflow PA. This agreement, along with the fact that the extension is only seen in 
the [SII] PSF of 2M1207A, is strong evidence against an artifact being the cause of the extension. 
Hence, we argue that we are in fact detecting the 2M1207A outflow in the [SII] image. However, 
note that no [SII] emission is detected at the source position in the UVES spectra (see Figure 
2). Why [SII] outflow emission should be detected in the FORSl image and not in the UVES 
spectra requires further investigation. Using Id flux calibrated spectra from the UVES pipeline, 
the flux of the 2M1207A [OI]/i6300 emission fine is estimated at 1 x 10"^^ erg/cm^s and the 



-8- 



continuum flux included in the range of the [SII] narrowband filter at 37 x 10 erg/cm^/s (0.6 x 
10~'^ erg/cmVs/A). Using UVES spectra of other BD/VLMS where both [01] and [SII] emission 
is detected we measured the average flux ratio of the [SII]i6716 / [OI]/16300 and [SII]/16731 / 
[OI]/16300 lines and use these values to make a realistic estimate of the ratio between the [SII] and 
[01] emission in 2M1207A. If this same ratio holds for 2M1207A, then the [SII] lines would have 
strengths of 0.18 and 0.3 x 10"'^ erg/cm^s, respectively. This puts the [SII] line at around the l-cr 
to 2-cr detection level, explaining its non-detection above the continuum. 

Due to the difference in critical density between the [OI]/16300 and [SII]/l/l6716, 6731 lines, 
the [01] emission is stronger than the [SII] emission at the source position and the [SII] emission is 
normally more extended. The critical densities of [OI]/16300 and [SII]/16731 lines, at T ~ 8,000 K, 



can be taken as 2 x 10^ cm~^ and 1.3 x 10"* cm~^ respectively ( Hartigan et al.||l995 1. Thus, knots 
further out in a jet are brighter in [SII] than [01] and this is illustrated, for example, for the Par- 
Lup3-4 outflow in Bacciotti et al. ( 2011| ). To check for faint extended [SII] emission we summed 



the spectra over +1" and the results are plotted in Figure 9. Indeed, once the spectra are summed 
in the spatial direction, both [SII] lines are detected at a velocity comparable to the [OI]/i6300 
velocity. The spectra have been smoothed in the same way as the [OI]/16300 spectrum and also 
shown is the Li I photospheric line which is used to measure the systemic velocity. Spectra from 
2006 and 2008 are compared and in both cases the [Sn]/I6716 line is brighter than the [SII]/16731 
line indicating a low value of n^ of ~ 1000 cm"^ ( Osterbrock||l988 1. The fact that some extended 



[SII] emission is detected in the UVES spectra supports our detection of the outflow in the image. 



3.3. [SII] Imaging: Extended Emission Features 

To the south-west of 2M1207A, a series of faint knot-like features which culminate in a bow- 
shock shaped object at a distance of ~ 16 " are detected (see Figures 7 and 10). These will now be 
referred to as features A to D, with D being the feature furthest from the BD. These sources seen 
in [Sn], R and I are in the direction of the red-shifted outflow and look very much "to the eye" 
hke they could be outflow knots. They are well fitted with a PA of 245° as expected for the red- 
shifted flow. However, it is possible that A - D could be faint stars or galaxies and therefore further 
analysis is needed to confirm that they are HH objects forming part of the 2M1207A outflow. 
Firstly a [SII]-Tband subtraction was carried out to isolate any outflow emission. In Figure 10 the 
FORSl [SII], R and I-band images and a SUSI2 I-band image downloaded from the ESO archive 
are compared. The features are clearly detected in all bands. The R-BESS filter covers the range 
5820-7320 A and the I-BESS filter the range 6990-8370 A. Hence, both filters include outflow 
tracers, especially the R-BESS which covers strong tracers like [OI]/16300, Hor, [NII]/16583 and 
[SII]/16731. Outflow tracers such as [FeII]/l7155 and [CaII]/l7291 are found in the range of the 
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I-band. However, as these lines are fainter than the lines in the range of the R-band and sometimes 
not present, the outflow emission in the I-band is generally significantly less than, for example, that 
of the [SII] narrowband filter ( Mundt et al.|1983 , Brugel et al.|1981| ). Therefore, I-band broadband 



images are routinely used to remove continuum emission from a narrowband image and thus to 



enhance any outflow emission. In Figure 11a Gaussian smoothed FORSl [SII]-I image of the 
region where it was possible to remove the I-band fringing is shown. The emission regions A-D 
are detected and the bow-shaped nature of feature D is obvious. The red line marks the position of 
outflow FA which well fits knot D in particular. The emission regions B and C are just above the 
noise in the [SII] -I image while A is the brightest. 

To further eliminate the possibility that A-D are stars / galaxies a [SII]-I versus R-I color- 
color analysis was carried out and is presented in Figure [12} A HH object is expected to display a 
blue broadband R-I color due to the contribution of the Halpha and [SII] 6716, 6731 lines to the 
R band, and an even bluer [SII]-I color caused by the larger fraction of the [SII] filter passband in 
which [SII] emission is noticeable, in contrast with the redder R-I and [SII]-I colors of continuum- 
dominated sources like stars and normal galaxies. Therefore this analysis should separate any HH 
objects from stars and galaxies. Firstly synthetic colors for stars with T^ff in the range 20000 K to 
1500 K are plotted (blue squares). These were calculated for the wavelength ranges of the FORSl 
filters using a simple blackbody model. The hottest stars lie close to the origin while the coolest 
BDs lie in the tail at the end of the sequence (upper right). The direction of the reddening is also 
plotted. In this figure the colors for both the TW Hya and Cha I star forming regions (based on the 
FOV of our images) are also plotted. The idea is to compare the colors of the sources in our images 
of the 2M1207 region not only to the synthetic colors but also to another star forming region and 
any known HH objects within that region. Cha I was chosen as suitable images were available in 
the ESO archive which included the protostar ESO-Ha 574 and its outflow ( jBacciotti et al.|20TT ). 



The black triangles represent the sources in the FOV of the TW Hya images and the black asterisks 
the sources in Cha I. For both star forming regions the sources follow the same sequence as the 
synthetic colors with the hottest objects clustering around the origin and the coolest in the tail. 

Looking first at the TW Hya data, the purple triangles are galaxies that could be identified in 
the field and the green triangles represent features A and D. B and C were omitted as they were too 
faint for their colors to be measured. While A lies amongst the stars and galaxies ruling out the 
possibility that it is a HH object, D is offset from the stellar sequence and is bluer than the stars 
and galaxies plotted. Thus considering the bluer nature of D, that it lies along the outflow V\ and 
that it is bow-shaped with the bow pointing away from the BD, we conclude that it is a HH object 
in the 2M1207A outflow. This conclusion is further strengthened by comparison with the knots 
in the ESO-Hor 574 outflow. ESO-Ha 574 powers a well-developed bipolar jet first observed by 



Comeron & Reipurth| ( |2006) . |Comer6n & Reipurth| ( |2006| ) identified for knots in the blue- shifted 



lobe of the ESO-Ha 574 flow which they also labelled A-D and Knot E in the red-shifted flow. 
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Here we label them HA574-A to HA574-D, to distinguish them from the features in the 2M1207A 
flow. In the images HA574-E was the only knot for which the colors could be measured as the 
others were not well enough separated and too faint in the I band. Its color is represented in Figure 
12 by the purple asterisk. HA574-E is also bluer and well separated from the stellar sequence as 
expected. The fact that it is bluer than D can be explained by the relatively weak [SII]/l/l6716, 
6731 emission from D compared to HA574-E. Flux calibrated spectra of the ESO-Ha 574 outflow 
which included HA574-A, HA574-B and HA574-E were also available (J. Alcala, priv. comm, 
also see Bacciotti et al. 201 1). These spectra were used to estimate the colors for these three knots 
(purple squares). The color for HA574-E lies close to that estimated from the images while there 
is a spread in the colors of HA574-A and HA574-B, again consistent with differences in shock 
conditions. 



4. Discussion and Conclusions 

Of the small number of BD outflows studied to date, 2M1207A is perhaps the most interesting 
due to its low mass, age of some 10 Myr and binarity. With the new results presented here we add 
considerably to our 2007 paper reporting the discovery of this outflow. The 2008 spectra confirm 
the original results, resolve for the first time for a BD the optical outflow, and allow an estimate of 
the outflow PA to be made. Furthermore, the outflow is detected in the [SII] image in the form of an 
elongation in the PSF and at least one HH object (feature D). Evidence in support of the extended 
features B-D being part of the 2M1207A outflow include: the fact that they lie along the estimated 
outflow PA; that the terminal feature D, is bow-shock shaped with the apex pointing away from 



the source and that D is separated in Figure 12 from the stars and galaxies. While the color-color 
analysis shown in Figure 12 demonstrates that at least D is a probable HH object, it shows that 
A is a star or galaxy. New higher S/N images would greatly add to current understanding of this 
outflow, as firstly one would be able to study the features with a much improved S/N and secondly 
a proper motion study, which would confirm the origin of B, C, could be conducted. Using the 



upper estimate of the disk inclination angle given by Riaz et al. ( 2012a|b ) (69°) and combining it 



with the measured radial velocities, we estimate the tangential velocity of the 2M1207A outflow at 
~ 20 kms"^ . While at the lower end of the range, this is in line with previous estimates of tangential 
velocities for a sample of BD outflows given by ,Whelan et al. ( |2009b ). Also the non-detection of 



a spectroastrometric signal for Ha is consistent with a low tangential velocity as with such low 
velocity one would expect only weak shocks and hardly any ionization of H in the jet and therefore 
no Ha. A velocity of ~ 20 kms"' would give feature D an age of some 173 years. Taking that 
any new images of 2M1207A would be observed in January 2013 at the earliest, this would mean 
a ~ 57 month time difference between these and the 2008 images. Thus a shift of at least ~ Qf.'S 
in features B-D would be expected, if they are indeed part of the outflow. Such proper motions 
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should be easily measured. The SUSI2 images shown in Figure 10 were obtained in January 2006, 
meaning that a shift of ~ 0'.'2, could reasonably be expected in the knots in the I band. Attempts 
to measure such a proper motion between these two sets of I-band images were inconclusive. A 
high quality spectrum along the PA of the outflow to search for shock related emission lines would 
also be of great benefit. New imaging and spectroscopic observations would also enable a search 
for emission features along the blue-shifted outflow to be conducted. Finally we note that while it 
is normally the case that the blue-shifted lobe is preferentially detected, there are cases where the 
red- shifted flow is more prominent (| Whelan et al.|2009b , Woitas et al.|2002 1. 



Results presented here also demonstrate that the 2M1207A outflow is well-coUimated. This 
confirms that BD outflows can be coUimated like low mass protostellar jets, further adding to 



the properties they share (Whelan et al. 2009b). As BD optical outflows had not been directly 
spatially resolved previous to this work, the extent of their coUimation was uncertain. However, 



the detection of a molecular outflow from ISO-Oph 102 by |Phan-Bao et al. (2008 1 did point to some 
degree of focussing, as it is a coUimated jet which is postulated to be driving the molecular outflows 
from low mass stars ( Bachiller||1996| ). Measurements of jet width as a function of distance along 
the jet for classical T Tauri stars (CTTS) have shown that coUimation occurs on scales of a few tens 
of AU from the star (Ray et al. 20071) which is consistent with magneto-hydrodynamical models 
( |Dougados|2009 ). As well as being coUimated, jets from low mass stars are also often episodic in 
that they are made up of a series of distinct knots assumed to be associated with separate accretion 
events. When discussing episodic outflows the striking H2 flow known as HH 212, driven by a 



Class young low mass star, immediately comes to mind (Zinnecker et al. 1998). The HH 212 
flow is very well defined with matching pairs of bow-shock shaped knots on either side of the 
driving source. By studying the individual knots in protostellar jets it is possible to get a picture 
of the mass loss and accretion history of the source. Our study strongly indicates that accretion 
and outflow activity proceeds in this extremely low mass object in an analogous way to low mass 
protostars, and by studying this flow we will be able to investigate the history of the 2M1207 
system. 

As the current sample of known BD outflows is small the present aim of our study of BD 
outflow activity is to greatly increase the sample. The approach of using SA to identify BD out- 
flows, while eff"ective, is slow as high quality, high S/N spectra are required and thus, observations 
with an 8 m class telescope or better are generally necessary. In addition, if one wants to thor- 
oughly investigate the properties of the outflow, including the morphology, kinematics and mass 
loss rate, spectra along the outflow PA are needed. Hence, to identify the outflow sources and 
constrain the outflow PA, two spectra at orthogonal PAs should be observed and then a follow-up 
spectrum must be taken along the estimated PA. We consider that the imaging observations of the 
2M1207A outflow presented here, ofl"er an underutilized method way of identify BD outflows. A 
[SII] imaging study to look for extensions in the PSFs of BD candidates would be a fast way of 
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identifying outflows as more than one source could be observed at the same time and we would 
expect to reach a better S/N in any [SII] outflow emission with an imager than a spectrograph for 
comparable exposure times. A serious limitation to imaging studies of jets from low mass proto- 
stars is the strength of the continuum emission which is many times stronger than the jet emission. 
For this reason the jet is often not seen in the image until several arc-seconds from the source. As 
the continuum emission in the vicinity of the [SII]/l/l6716, 6731 lines in BDs and VLMSs is fainter 
than in CTTSs, this enhances the possibility of detecting any outflows in an imaging study. 

Furthermore, a significant number of the BDs / VLMSs shown spectroscopically to have 
outflows (4 out of 6) also have near-edge-on disks CBacc iotti et al.||2011 ; |Whelan et al.||2009b|a ) . 



The eff'ect of this is to obscure the emission from the BD or star and further enhance the contrast 
between the jet and continuum emission in favor of the jet emission. For example, the 100 Mjup 
mass object Par-Lup3-4 is known to be under-luminous and an examination of archived FORS 1 
images and UVES spectra allows the ratio of the [SII] jet emission to continuum emission in the 
range of the [SII] -1-62 filter to be estimated at 50%. Thus, the Par-Lup3-4 jet is easily detected 



in the FORS image (Fernandez & Comeron 2005). We are confident of repeating the imaging 



observations of 2M1207A discussed here for other BDs and especially those with edge-on disks. 
Such observations would of course immediately give the jet PA allowing for a more detailed follow- 
up study of accretion and outflow activity from high quality spectra. 
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Fig. 1.— FORSl R-band image of 2M1207 with the UVES slits and sht PAs shown. The shts 
are drawn to scale. The spectro-astrometric offsets measured at the two slit R^s, for the blue and 
red-shifted [OI]/16300 emission, are represented by the blue and red arrows. The size of the arrows 
do not represent the magnitude of the measured offsets. However, they accurately represent the 
relative magnitudes between the two observations, in that the offsets measured at R^ = 80° are 
approximately twice those measured at PA = 0°. 
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Fig. 2. — The full spectral range of the UVES observations aimed at studying the 2M1207A 
outflow. This ID spectrum was extracted from the source position using the ESO UVES pipeline 
and was taken on May 16th 2006. Hor and HjS are by far the strongest lines and while [OI]/i6300 
Une is detected no [SII] emission is detected above the noise at the source position. In the inset 
showing the region of the [SII]/l/l6716, 6731 lines, the dashed red lines mark the kms"^ position 
for the two emission regions. 
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Fig. 3.— Left: Normalized [OI]i6300 line profile from the 2008 (solid line) and 2006 (dashed line) 
spectra. The line region has been smoothed to increase the S/N, as outlined injWhelan et al. (2007 1 
and in both cases the line is double peaked. Right: Position- Velocity diagram of the [OI]/16300 
line region from the 2008 spectra. The line region is spatially resolved and clearly originates in a 
bipolar outflow. The contours begin at 3 times the r.m.s noise and increase in intervals of the r.m.s 
noise. The zero spatial offset line is measured from an accurate mapping of the continuum position 
using SA. The [OI]/16300 night sky line has been subtracted from the spectrum. 
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Fig. 4.— Left: SA of the 2008 [OI]/16300 line. The SA clearly confirms what is revealed in the PV 
diagram in Figure [3} that the blue and red peaks are offset in opposing directions. The dashed line 
marks the l-cr error in the centroid fitting. Right: The spectro-astrometric signatures measured in 
the [OI]/16300 line in 2006 and 2008 are combined to constrain the outflow PA. Note that the x and 
y axes do not represent orthogonal slit PAs but rather PAs of 0° and 80°. 
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Fig. 5. — Spectro-astrometric analysis of the Ha emission in all the UVES spectra obtained for the 
outflow study. The purpose here was to check if an outflow component could be detected as the 
Ha line varied. No off'set is detected in any of the spectra hence it can be concluded that the vast 
majority of the Ha emission is likely tracing accretion in the BD and that the Ha line can be used 
to provide a reasonable estimate of the mass accretion rate. Note that the line is highly variable 



as outlined by Scholz et al. (2005). The appearance of a small bump between the primary and 
secondary peaks was not previously detected in any variability study. 
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Fig. 6. — Comparing the Hfi line profiles in the 2006 spectra. Both the line and continuum vary. 
No spectro-astrometric signal was detected for the Hfi line. 
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Fig. 7. — [SII] image in the region of 2M1207. A zoom on 2M1207 (the A and B components to 
the system are not resolved) is shown in an inset. Note the slight elongation in the PSF along a 
similar PA to the outflow PA derived from the spectro-astrometric analysis. The PA of 2M1207B 
with respect to 2M1207A is 125°. The position of the possible outflow features are marked. As 
described in the text two images with exposure times of 1800s were obtained. The image with the 
better seeing is shown here. 
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Fig. 8.— Comparison of the PSF of 21M1207A in [SII] with the PSFs of 2M1207A in R and the 
PSFs of two nearby stars in [SII]. The purpose here is to demonstrate that the 21VI1207A PSF in 
[SII] is elongated along the outflow PA with respect to the other PSFs. An elliptical fit to the PSFs 
is shown in the inset. The ellipticity of the fits shows that the 21VI1207 [SII] PSF is indeed more 
elliptical and has a PA estimated from the fit at 62° + 6° . This PA is compatible with the outflow 
PA as estimated using SA. We use this comparison as strong evidence that we are detecting the 
2M1207A outflow in the [SII] image. 
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Fig. 9.— The Li I absorption line and the [SII]/1/16716, 6731 doublet in the 2006 (left column) and 
2008 (right column) spectra. The Li I spectrum is extracted from the source position and used to 
estimate the systemic velocity at ~ +2 kms ^ The [SII] lines are not detected at the source position 
but only by summing over + 1 ". As explained in the text this is likely due to the lower critical 
density of the [SII] lines compared to the [OI]/16300 line. 
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Fig. 10. — Here the knot-like features seen to the S-W of 2M1207A with a PA of ~ 245° are 
compared in the [SII], R-band and I-band images. Two of the emission features are detected in 
the FORS 1 I-band image. Also shown in a SUSI2 image taken two years previous to the FORS 1 
observations 
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Fig. 11. — Gaussian smoothed [SII]-I image. 4 separate emission features are revealed with the 
final one having a clear bow-shock shape. The positions of the features are well fitted with a PA of 
245° which is the same PA as estimated from the spectro-astrometry. 
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Fig. 12. — Using a color-color analysis to separate HH objects from stars and galaxies. The aim 
of this analysis is to compare the colors of the objects in the TW Hya images with synthetic colors 
and with an analysis of the protostar ESO-HA 574 and its associated HH objects, in order to test 
if the faint knot-like features A and D are also HH objects. The squares represent synthetic colors, 
the circles colors estimated from spectra, the triangles represent objects in TW Hya included in 
the FOV of the 2M1207 images and the asterisks represent objects in Cha I included in the FOV 
of the ESO-HA 574 images. The red arrow is the reddening vector. Specifically the blue squares 
are synthetic colors for blackbodies in the range 20000 K to 1500 K, calculated in the range of 
the FORSl filters. The purple circles are the color estimates for the knots HA574-A, HA574-B 
and HA574-E, calculated from XSHOOTER spectra. The black asterisks are stars in Cha I and 
the purple asterisk marks the color of HA574-E measured from the image. The black triangles 
are stars in the TW Hya region, the purple triangles are identified galaxies and the green triangles 
features A and D. The marked errors in the colors for HA574-E and A, D are quite large while the 
errors for the other points are found to be < 0.02. The position of feature D, the bow-shock shaped 
feature lies with the other known HH objects while A lies amongst the stars and galaxies. Based on 
this analysis and other evidence discussed we conclude that Z) is a shock in the 2M1207A outflow. 



